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I . INTRODUCTION 
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Although  a great  deal  of  progress  has  been  made  in  describing 

the  gaseous  and  solid  states  of  matter,  the  situation  with  liquids  is 

(a) 

far  from  clear  and  highly  disputed  [1*2].  ' Of  all  the  properties  of 

liquids,  viscosity  has  probably  received  the  greatest  attention  [67]. 
The  viscosity  data  that  have  been  collected  on  mixtures,  however, 
usually  have  not  covered  both  full  composition  ranges  and  a range  of 
temperatures.  The  present  work  provides  such  data  for  a variety  of 
systems  that  range  from  ideal  solutions  in  the  thermodynamic  sense,  to 
ones  in  which  two  liquid  phases  are  formed.  Viscosity  data  at  temper- 
atures and  compositions  in  the  neighborhood  of  phase  separation  have 
been  particularly  lacking. 

With  the  large  amounts  of  data  available,  numerous  attempts 
have  been  made  to  develop  theories  of  viscosity.  These  theories  have 
been  collected  by  Partington  [71]  and,  more  recently,  in  a thorough 
review  by  Brush  [lf>].  In  many  of  these  cases  the  theories  have  been 
developed  using  the  data  for  only  a few  compounds.  This  suggested 
that  a further  examination  should  be  made  to  determine  which  of  these 
theories  are  simply  empirical  fits  of  data  and  which  offer  a reasonable 
description  of  the  liquid  state  on  the  molecular  level. 

(a) 

v 'Numbers  in  square  brackets  [ ] refer  to  references  in  the 
Bibliography. 
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In  the  present  work  some  of  the  factors  involved  in  flow  have 
been  separated.  These  include  a volume  factor  and  a local  inhomogen- 
iety  factor.  It  is  also  proposed  that  the  interaction  between  rigid 
acentric  molecules  affects  the  molecular  rotational  motion.  This  effect 
is  accompanied  by  a rotational  energy  of  activation  for  viscous  flow. 
Ibis  may  explain  many  of  the  viscosities  and  the  temperature  coeffi- 
cients of  viscosity  which  deviate  widely  from  current  theories. 

Some  preliminary  work  on  surface  tensions  has  been  done  and 
a unique  observation,  that  of  liquids  rising  in  a thin  film  as  much  as 
8 cm.  on  a glass  wall,  has  been  noted  for  certain  hydrocarbon- 
fluorocarbon  mixtures. 


II . APPARATUS 


Density.  Pycnometers  described  by  Lipkin  [57]  were  calibrated 
with  distilled  mercury  at  each  temperature  for  which  liquid  densities 
were  determined.  The  nominal  capacity  of  the  pycnometers  was  1 milli- 
liter and  10  milliliters.  The  precision  of  the  density  values  obtained 
is  ± 0.0002  grams  per  milliliter  in  the  one-milliliter  pycnometers. 

Viscosity . The  Cannon-Fenske  modification  [5]  of  the  Ostwald 
viscometer  was  calibrated  with  an  oil  of  standard  viscosity  obtained 
from  the  National  Bureau  of  Standards  and  with  a sample  of  pure  toluene. 
The  physical  properties  of  the  toluene  (density,  boiling  point,  and 
refractive  index)  were  identical  with  those  given  in  "Selected  Values 
of  Physical  and  Thermodynamic  Properties  of  Hydrocarbons,"  published  by 
Project  1|1|  of  the  American  Petroleum  Institute  at  the  Carnegie  Institute 
of  Technology.  The  viscometers  were  calibrated  at  each  temperature  for 
which  viscosities  were  determined.  Ubbelohde  suspended  level  viscom- 
eters [5]  calibrated  by  the  Cannon  Instrument  Company  were  used  for 
about  60  per  cent  of  the  measurements  and  reproduced  the  data  obtained 
with  the  other  viscometers.  The  clocks  were  calibrated  with  reference 
to  the  standard  time  intervals  transmitted  by  radio  station  WWV.  The 
precision  of  the  viscosities  is  2 parts  in  1000. 

Temperature . All  mercury-in- glass  thermometers  and  all  thermo- 
couples used  in  measuring  thermostat  temperatures  and  boiling  tempera- 
tures (except  those  reported  in  distillations)  were  calibrated  against 
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an  iron-constantan  thermocouple  calibrated  at  the  National  Bureau  of 
Standards.  The  accuracy  of  the  standard  thermocouple  was  ± 0.01°C. 

Vapor  Pressure.  The  vapor  pressure- temperature  determinations 
on  the  pure  materials  were  obtained  by  boiling  the  liquid  in  a Cottrel 
apparatus  under  a measured  air  pressure.  The  boiling  temperatures  were 
Ouv/aj.ned  xrom  calibrated  thermocouples.  The  pressure  was  measured  with 
a precision  of  ± 0.05  mm.  Kg  with  a one-meter  cathotometer  and  a mer- 
cury column  20  mm.  in  diameter. 

Refractive  Index.  Abbe- type  instruments  were  used  for  these 
measurements.  One  refractometer  had  a range  of  1.2  to  1.1*  while  the 
other  had  the  standard  range  of  1.3  to  1.6. 

Vapor  Chromatograms.  A Perkin-Elmer  Model  !$k  "Vapor  Fractom- 
eter"  was  used  for  the  chromatographic  determinations  of  purity.  The 
Ccir'r^-er  §ases  used  were  nitrogen  and  helium.  The  stationary  media  used 
to  separate  fluorocarbon  materials  were  those  reported  [77]  to  be  effec- 
tive (e.g.,  Cl(CF2CFCl)3CF2COOC2H^  and  (CjF  ) N). 

Sar^ace  tensions.  The  surface  tensions  were  measured  by  a drop 
weight  technique  [93].  The  apparatus  consisted  of  a glass  capillary 
tube  with  a tip  which  was  machined  to  a circular  diameter  of  2.1*6  ± 0.01 
mm.  and  polished  flat.  Drops  were  slowly  formed  on  this  tip  (one  to 
four  minutes  per  drop]  and  allowed  to  fall  into  a small  weighing  bottle. 
The  weighing  bottle  was  located  upon  a small  platform  in  an  enclosure 
containing  a sufficient  amount  of  the  liquid  being  tested  to  insure 
vapor  saturation  of  the  atmosphere. 

Between  ten  and  twenty  drops  were  used  in  each  determination 
and  each  ol  the  determinations  was  duplicated.  An  average  precision 
Ox  approximately  1 per  cent  was  obtained. 
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An  apparatus  has  been  designed  which  will  give  unusually  high 
accuracy  in  measuring  relative  surface  tensions  although  no  measurements 
with  ohis  equipment  are  included  in  this  work.  This  apparatus  consists 
of  a U tube  in  which  one  arm  is  a capillary  and  the  other  is  a large 
precision  bore  glass  tube.  The  large  tube  is  equipped  with  a ground- 
glass  cover  to  prevent  evaporation  of  the  liquid.  Liquid  of  known  den- 
sity is  admitted  to  the  large  tube  until  the  liquid  level  in  the  capil- 
lary reaches  an  etched  mark.  The  difference  in  weight  of  the  apparatus 
before  and  after  filling  allows  the  capillary  rise  to  be  calibrated  very 
accurately.  Since  each  liquid  is  added  until  the  level  in  the  capillary 
has  reached  the  same  point,  the  effect  of  any  irregularities  in  the  bore 
or  optical  aberrations  in  the  glass  are  eliminated. 


Ill . MATERIALS 


Perfluoroheptane,  C.,F,^.  A fluorocarbon  sample  (Minnesota 
Mining  and  Manufacturing  Company  material  FM-3130)  was  shown  by  vapor 
chromatography  to  contain  at  least  19  components,  the  major  component 
(approx.  73%)  being  C„F.^.  A charge  of  1319  grams  of  this  material  was 

distilled  in  a column  equal  to  60  theoretical  plates  when  distilling 

(a) 

hydrocarbons.  The  constant  boiling  center  portion  of  the  product, 
639-2  grams,  gave  only  two  peaks  upon  analysis  by  vapor  chromatography. 
The  second  peak  constituted  approximately  1 to  2 per  cent  of  this  mater- 
ial. This  constant  boiling  material  was  then  azeotroped  with  198.7 
grams  of  n-heptane  in  the  60-plate  column.  An  improved  material  came 
overhead  as  the  initial  product  in  the  form  of  an  azeotrope  boiling 
approximately  9C°  below  the  flat  obtained  in  the  original  distillation. 
The  first. 387 -3  grams  of  this  azeotrope  were  cooled  to  dry-ice  tempera- 
ture where  the  hydrocarbon  and  fluorocarbon  are  almost  completely  immis- 
cible. The  two  layers  were  separated  mechanically.  This  gave  336. b 
grams  of  fluorocarbon- rich  layer  which  was  placed  in  a distillation 
column  equal  to  30  theoretical  plates  when  distilling  hydrocarbons. 

/ \ 

\ di  J 

yThe  distillation  columns  used  in  this  work  were  packed  with 
1/16  inch  single- turn  helices.  It  has  been  shown  [97]  that  such  col- 
umns have  about  one-half  as  many  theoretical  plates  when  distilling 
fluorocarbons  as  when  distilling  hydrocarbons.  The  values  quoted  are 
for  the  number  of  theoretical  plates  in  the  columns  when  used  for  dis- 
tilling hydrocarbons. 
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ihc  remaining  traces  of  dissolved  n-heptane  were  then  removed  as  the 

azeotrope.  When  analysis  by  vapor  chromatography  showed  no  more 

n-heptane  in  the  overhead  the  distillation  wTas  stopped.  The  residue 

from  this  distillation  (322.9  g.)  was  taken  as  the  pure  C Jib  . for  the 

7 16 

viscosity  studies. 

under  exoremely  high  sensitivity  on  the  vapor  chromatograph, 

V helium  carrier  gas)  this  C^F^  material  showed  a small  shoulder  on  the 

back  of  the  main  peak.  Since  these  two  materials  appear  so  close 

together  it  has  been  assumed  that  they  are  two  different  C„Fn.  isomers. 

7 16 

m_s  i s further  borne  out  Dy  the  nuclear  magnetic  resonance  spectrum 

run  on  another  C^F^  sample  from  the  same  source  which  gave  an  almost 

ice.. ^j_cal  trace  on  tne  vapor  chromatograph  under  the  same  high  resolution. 

The  second  component  showed  the  presence  of  (CF^)2  C structure 

in  an  amount  no  greater  than  10  per  cent  as  a C-Fl6  isomer.  Askew  [U] 

has  identified  this  isomer  as  2-CF^C^F^.  The  only  other  peaks  in  the 

nuclear  magnetic  resonance  (n.m.r.)  spectrum  were  those  for  n-C  F 

7 16 

This  C7Fl6  material  has  a boiling  point  of  82.3°C  at  760  mm.  Hg. 
Perfluorocyclic  Ether,  CgF^O.  Fractionation  of  29^1  grams  of 
Minnesota  Mining  and  Manufacturing  Company  Fluoro chemical  101  in  a 60- 
pxaue  distillation  column  yielded  617.3  g.  of  material  refluxing  at 
101° C with  a refractive  index  of  from  n2^  = 1.2772  to  n2^  = 1.2771;. 

\ apor  chromatography  showed  this  material  to  contain  one  major  component 
(93$)  and  three  minor  components.  This  material  was  then  refluxed  77 
hours  with  80  g.  of  K0H  to  decompose  any  hydrogen-containing  material. 

At  the  end  of  this  time  the  K0H  pellets  had  turned  dark  brown.  The 
fluorocarbon  was  removed  and  again  refluxed  with  80  g.  of  fresh  K0H 
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for  63  hours.  At  the  end  of  this  treatment  the  fluorocarbon  was  washed 
four  times  with  equal  volumes  of  distilled  water,  dried  with  CaC^,  and 
simple  distilled.  This  product  had  a refractive  index  of  = 1.2775. 
Vapor  chromatography  still  indicated  the  presence  of  small  amounts  of 
three  minor  components.  The  568  grams  of  material  that  remained  from 
the  above  treatment  was  distilled  in  the  60-plate  column  with  127  g. 
of  n-heptane.  Most  of  the  remaining  impurities  came  overhead  in  a frac- 
tion refluxing  at  82.1  to  83.9  C.  These  latter  fractions  were  combined 
with  the  residue  and  cooled  to  dry-ice  temperature.  The  bottom  layer 
contained  k2 9 g.  This  layer  was  distilled  in  a 30-plate  column  to 
remove  the  final  traces  of  n-heptane  as  the  azeotrope.  The  residue 
showed  a small  shoulder  on  the  back  of  the  main  peak  under  high  sensi- 
tivity in  the  vapor  chromatograph.  This  impurity  represents  less  than 

1 per  cent  of  the  sample.  This  material  was  taken  as  the  CoF-,.0  for  the 

o 16 

viscosity  studies.  It  has  a boiling  point  of  103. 0°C  at  760  mm.  Hg. 

C8F16°  a cyclic  perfluoro  ether  of  unknown  structure.  It  is 
believed  that  the  oxygen  atom  is  in  the  ring,  however  [66]. 

2,2, 3-Trichloroheptaf luorobutane , C,  Cl^  and  1, 2-Dichlorohexa- 
fluorocyclopanter.e-1 , C^Cl^  . These  materials  used  in  the  viscosity 
studies  were  constant  boiling  fractions  obtained  from  the  distillations, 
in  a 60-plate  column,  of  two  crude  materials  obtained  from  the  Hooker 
Electrochemical  Company.  Vapor  chromatography  indicated  only  a very 
small  amount  of  impurity  (less  than  0.5%)  in  both  materials. 

n -Hexane,  n-C^H^.  Commercial  grade  hexane  was  washed  with 

H2S0, , passed  through  silica  gel,  and  distilled  in  a 60-plate  column. 

±he  inactions  containing  predominantly  n-hexane  were  again  washed  with 
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anc^  aze°Vroped  with  excess  ethanol  in  the  60-plate  column.  The 
azeotrope  was  washed  with  water,  dried,  and  the  remaining  ethanol 
removed  as  the  azeotrope  in  a 30-plate  column.  The  n-hexane  used  in 
these  studies  was  a constant  boiling  fraction  from  this  distillation 
in  which  no  impurities  were  indicated  by  vapor  chromatography. 

b00>-  i-'&n5 , 2 , d. , I4-1  rimethylpentane , i-CgH^g.  The  isooctane 
used  m these  studies  was  Phillips  99  mole-per  cent-minimum  grade. 
Vapor  chromatograms  of  this  material  showed  no  impurities. 

The  physical  properties  of  these  compounds  are  summarized 

in  Table  1. 


Physical  Properties  of  Compounds  Used  in  Experimental  Work 
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TABLE  1 (Continued) 
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IV.  EXPERIMENTAL  RESULTS 


Viscosity  Data.  The  viscosities  have  been  studied  for  ten 
binary  systems.  Viscosities  and  densities  were  measured  for  eleven 
compositions  in  each  system  at  3C°  intervals  from  2£°C  to  li£°C.  The 
experimental  kinematic  viscosities  are  presented  in  Figures  1 through 
10.  These  figures  are  all  drawn  to  the  same  scale.  The  composition  in 
each  mixture  is  expressed  in  terms  of  the  volume  fraction  at  25°C  of  the 
component  with  the  greater  kinematic  viscosity  as  a pure  compound. 

This  series  of  systems  allows  a systematic  study  to  be  made  of 
soxutions  varying  from  the  ideal  to  those  with  large  positive  deviations 
in  ^ne  activity  coefficient  from  the  ideal  solution  behavior.  These  sys- 
tems are  listed  in  Table  2 in  the  order  of  increasing  positive  deviation 
from  thermodynamically  ideal  behavior.  The  K values  defined  in  this 
table  are  those  which  appear  in  Hildebrand's  theory  [1|1].  In  this  theory, 

tne  logarithm  of  the  activity  coefficient  of  compound  "a"  is  K cp^/T 

aTb  3 

where  cpb  is  the  volume  fraction  of  component  b.  K equals  zero  for  an 
ideal  solution  and  increases  with  larger  deviations  from  the  thermo- 
dynamic xdeal  solution  behavior . Solutions  of  fluorocarbons  with  hydro- 
carbons require  special  treatment  [76,  98].  The  K values  listed  in 
Table  2 have  been  shown  [98]  to  describe  the  activity  coefficients  of 
systems  one  through  four  fairly  accurately.  From  Tables  3 through  12 
it  can  be  seen  that  the  volume  change  on  mixing  increases  rather  steadily 
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Figure  1. — Viscosities  of  the  C„F-, , 

I 16 


C8Fi60  System 


Ill 


Figure  2. -Viscosities  of  the  C^Cl^  - CjCl^  System 
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Figure  3. — Viscosities  of  the  C^Cl^F^  - CgF-^0  System 
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Figure  i;. — Viscosities  of  the  C^Cl^F^  - CgF.^0  System 
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Figure  — Viscosities  of  the  C^Cl^F^  - C^F^  System 


Viscosity  in  Centistokes 
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Figure  6. — Viscosities  of  the  n-C^H^  - C^Cl^F 


6 System 


Viscosity  in  Centistokes 


Figure  7. — Viscosities  of  the  C^Cl^F^  - i-CgH^g  System 


Viscosity  in  Centistokes 
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Figure  9. — Viscosities  of  the 


--C6HlU 


CgFisO  System 
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Figure  10.— Viscosities  of  the  i-CgH^  - System 
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TABLE  2 

Systems  Used  in  Experimental  Work 
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solubility  parameter  [1*1]. 


2h 


from  a very  small  negative  value  on  the  order  of  -0.01  per  cent  for 

system  one,  through  larger  positive  values  up  to  approximately  1 to 

<# 

3 per  cent  of  the  total  volume  for  system  ten. 

Vapor-liquid  equilibrium  data  for  systems  one  and  two  show  that 
these  components  form  solutions  very  close  to  ideal  while  system  three 
is  sufficiently  non-ideal  to  form  an  azeotrope  at  atmospheric  pressure 
[98].  Vapor-liquid  equilibrium  data  are  not  available  for  system  five, 
but  for  system  four  positive  deviations  from  Raoult's  law  are  known. 
Systems  five,  six  and  seven  are  shown  to  be  less  ideal  than  those  of 
the  chlorofluorocarbons  with  fluorocarbons  by  the  K- values  of  Table  2. 
(It  is  known  that  C^Cl^F^  forms  a minimum  boiling  azeotrope  with 
n-heptane  at  atmospheric  pressure  [98].)  No  thermodynamic  data  are 
available  for  system  eight,  but  like  all  fluorocarbon-hydrocarbon  solu- 
tions this  one  undoubtedly  has  large  activity  coefficients. 

No  thermodynamic  data  are  available  for  system  nine.  The  crit- 
ical unmixing  temperature  was  determined  in  this  work  at  26°C  at  27 
mole  per  cent  CgF^O. 

For  system  ten,  both  vapor-liquid  [69]  and  liquid-liquid  [i|0] 
data  are  available.  The  kinematic  viscosity  of  this  system  (Figure  10) 
has  been  examined  in  greater  detail  than  that  of  the  above  systems. 

As  the  temperature  is  decreased  from  ii5°C  to  25>°C,  the  portion  of  the 
curve  in  the  middle  composition  range  becomes  at  first  linear  at  35° C, 
and  then  exhibits  a slight  rise  at  30°C  which  finally  becomes  a max- 
imum at  25°C.  The  critical  unmixing  temperature  for  this  system  has 
been  observed  variously  at  23.7°C  in  reference  [1|0]  (for  n-C^F^)  and 
at  23.65°C  in  this  work  for  the  C-F-^  material  used.  The  dashed  line 
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Assuming  no  volume  change  on  mixing. 


TABLE  11 

Density , Volume  Change  on  Mixing  and  Refractive  Indices 
for  n-C^H^  - CgF^-0  System 
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33 


0) 

H O 

CA 

vO 

o 

CM 

CA 

o O 

CM 

C\i 

o 

On 

<A 

CO 

e m 

• 

• 

• 

• 

• 

\ -4 

CA 

CA 

pj 

fA 

CA 

r — 1 

• 

■d 

+ 

bO 

c 

•3  o 

vO 

\T\ 

3 

oo 

CO 

CM 

e o 

CO 

vO 

O 

r- 

1A 

UA 

• 

• 

• 

4 

• 

• 

> CM 

<3 

CM 

+ 

C\J 

<A 

<A 

CM 

rH 

fA  O 

o o 

rt  h 

c^- 

co 

CO 

UA 

O 

O 

ca 

CM 

1A 

CA 

rH 

e 

• 

• 

• 

• 

• 

• 

1 I73  bO 

C\J 

C\J 

CM 

rH 

rH 

i — 1 

o 

o 

*UA 

CM 

-P  -« 

o 

*LA 

1A 

fA 

vO 

CA 

•»  *rl  • 

rH 

CM 

CO 

Pt 

On 

O 

o.  CO  H 

CA 

CM 

O 

vO 

"LA 

vO 

C!  S 

vO 

*LA 

Ov 

CO 

C"~- 

<D  \ 

• 

• 

• 

CO 

• 

• 

• 

O 

rH 

rH 

rH 

<1) 

o 

o 

o 

bO 

CO 

cti 

pu. 

TJ 

•H 

. 

1 

o *n 

rA 

0\ 

UA 

■H 

CM 

o 

o 

rt  X O 

CM 

CO 

M3 

PI 

1A 

vO 

OO 

U Q)  O 

CO 

CO 

Ov 

<A 

Pt 

1A 

Ch  "O  1A 

C\J 

CM 

CM 

o 

CA 

<A 

fA 

(0  q cm 
CtJ  *H 

• 

• 

• 

& 

• 

• 

• 

rH 

rH 

rH 

E-i 

rH 

rH 

rH 

pf 

rH 

w 

VO 

o 

rH 

CA 

rH 

Ov 

CM 

CA 

vO 

On 

rH 

a)  1. 

ON 

Pt 

Pt 

fA 

O 

Ov 

C^- 

o 

rH  q | 

• 

• 

• 

• 

• 

• 

• 

o 

vO 

O 

fA 

Pt 

pf 

CA 

o 

C^- 

Pt 

s 'ba 

rH 

CA 

Pt 

UA 

vO 

r- 

CO 

co 

On 

o -nt 
rt  i — i 

v-'  W O 

<0  M3  O 

6 0 UA 

» 

O 

vO 

fA 

"UA 

VO 

1A 

Pt 

3 1 CM 

O 

Ov 

O 

1 — f 

rH 

<A 

o 

f- 

rH  C| 

• 

• 

• 

• 

• 

• 

. 

• 

• 

O -P 

o 

o 

ON 

O 

o 

o 

O 

o 

On 

> ba  rt 

rH 

CM 

CM 

m 

vO 

a- 

CO 

. CO 

Assuming  no  volume  change  on  mixing. 


TABLE  12 

Density,  Volume  Change  on  Mixing  and  Refractive  Indices 
for  i-CgH^g  - System 
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on  Figure  10  is  a portion  of  the  approximate  locus  of  the  viscosities 
corresponding  to  the  consolute  curve.  Me  is  the  viscosity  at  the  crit- 
ical temperature  of  23.65°C. 

The  isotherms  of  viscosity  versus  composition  for  system  ten 
(Figure  10)  are  similar  to  those  of  system  nine  near  the  unmixing  curve. 
The  viscosities  of  systems  eight,  nine,  and  ten  reveal  that  changes  in 
the  structure  of  the  solution  which  ultimately  produce  separation  into 
two  equilibrium  phases  begin  at  temperatures  at  least  10C°  above  the 
temperature  at  which  two  liquid  phases  are  observed. 

Surface  Tension  Data.  Surface  tensions  were  determined  for 
eight  of  the  binary  systems  (all  except  systems  four  and  five).  The 
values  for  each  system  were  obtained  at  three  different  temperatures 
(25°C,  35°C  and  U^°C ) and  at  six  compositions.  These  data  are  presented 
in  Figures  11  through  18. 

The  deviations  from  linearity  in  volume  fraction  is  seen  to 
parallel  closely  the  deviations  from  linearity  of  the  kinematic  vis- 
cosities for  the  same  systems.  Although  these  surface  tension  values 
are  not  extremely  accurate  (deviations  up  to  3 per  cent  from  reported 
[1]  values)  the  departure  from  linearity  is  well  established.  No  further 
attempts  were  made  to  correlate  the  surface  tensions  with  the  viscosities 
when  it  was  observed  that  neither  the  magnitude  nor  the  temperature  coef- 
ficient of  surface  tension  paralleled  the  corresponding  viscosity  values 

♦ 

for  the  pure  components.  It  should  be  noted,  however,  that  any  adequate 
theory  of  the  liquid  mixtures  should  explain  both  surface  tension  and 
viscosity. 


J 


36 


J 


Figure  11.— Surface  Tensions  of  the  C^F.^  - CgF^O  System 
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Volume  Per  Cent  C^Cl^F^  at  25°C 
Figure  13.—  Surface  Tensions  of  the  C,  Cl  F - CfiF  .0  System' 
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Figure  li*. — Surface  Tensions  of  the  n-C^H^  - C^Cl^F^  System 
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Figure  15.— Surface  Tensions  of  the  C^Cl^  - i-CgHl8  System 
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Figure  16. — Surface  Tensions  of  the  i-CgH^g  - CgF^O  System 
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Volume  Per  Cent  CgF^O  at  25°C 
Figure  17.—  Surface  Tensions  of  the  n-C^^  - CgF.^0  System 
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Relative  Surface  Tensions  by  a Differential  Capillary  Rise 
Method.  Due  to  the  lack  of  both  precision  and  accuracy  encountered  in 


the  drop-weight  surface  tension  determinations,  a new  apparatus  has  been 
designed  which  will  permit  the  measurement  of  relative  surface  tensions 
with  a high  degree  of  accuracy. 

The  major  problem  in  the  drop-weight  method  is  that  the  size  of 
the  drop  is  dependent  upon  its  rate  of  formation.  Similar  problems  are 
present  in  any  procedure  in  which  the  surface  area  is  varying  at  the 
time  of  measurement.  While  the  capillary  rise  method  is  free  from  this 
objection,  several  other  difficulties  are  encountered.  In  the  apparatus 
designed  in  this  work  three  of  the  major  sources  of  error  are  eliminated. 

Jones  and  Ray  [5l]  have  developed  an  equation  for  surface  ten- 
sions using  the  capillary  rise  method: 


o (D  -P  ) 

c _ c c' 

a (D  -P  ) 

o v o o' 


1 + 


1 

rrR2ho(l-r/R) 


where : 


c = surface  tension 
D = density  of  liquid 

0 = density  of  gas  phase  (air  and  vapor) 

W = weight 

R = radius  of  large  tube 

r = radius  of  capillary  tube 

h = true  capillary  rise  above  infinite  surface 


subscripts: 

o = reference  liquid 
c = test  liquid. 


The  above  equation  does  not  assumes 


1.  Infinite  reservoir 

2.  Any  particular  shape  of  meniscus 

3.  Uniform  capillary 

is.  Optically  uniform  apparatus. 

Jones  and  Ray  measured  hQ,  r and  R independently.  This  can, 
however,  be  avoided  by  writing  the  equation  as 


— = constant  (at  each  T) 

t®  h (1-r/fo) 


fa  (D  -3  ) 

c o o 

a (b  -p  ) “ 

L ov  C c' 


This  constant  may  be  evaluated  if  two  standards  are  available  (substi- 
tute for  a^  and  a^  in  the  last  equation). 

In  practice  enough  material  is  added  to  the  large  tube  to  bring 
the  liquid  level  in  the  capillary  exactly  to  an  etched  mark.  This  amount 
of  material  is  determined  by  the  weight. 

Do,  3q  and  Wq  are  constants  which  must  be  evaluated  at  each  T. 

aQ  is  a known  constant  or  ^ may  be  used  to  give  the  relative 

o 

surface  tensions. 


D , 3 and  W must  be  determined  for  each  surface  tension  exper- 
c c c 


iment . 


R must  be  of  uniform  inside  bore  over  the  working  range.  Such 
tubes  of  extremely  uniform  bore  are  commercially  available. 

The  liquid  height  in  the  capillary  must  be  observed  from  the 
same  height  and  angle  so  that  any  optical  distortions  are  compensating 
and  the  apparatus  must  be  aligned  in  the  bath  in  the  same  manner  for 
each  determination.  A permanently  mounted  support  can  assure  this 
alignment . 
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Liquid  Rise  on  Glass  Walls.  One  of  the  most  unusual  properties 
of  the  systems  studied  is  that  certain  of  the  mixtures  will  ''creep"  as 
a thin  film  up  the  glass  walls  of  the  apparatus  in  a manner  highly  sug- 
gestive of  the  well-known  super- fluidity  of  liquid  helium. 

Of  the  binary  mixtures  studied,  systems  seven  (C^Cl^  F6  - i-C8Hl8>  > 

eight  (i-CgHl8  - > and  ten  (C^F^  - i-CgH^g)  exhibited  this 

unusual  behavior.  This  phenomenon  was  observed  in  the  hydrocarbon- rich 
mixtures  of  these  systems.  Isooctane  is  the  only  material  common  to 
all  three  of  these  systems,  and  as  a pure  material,  it  does  not  creep 
up  glass  walls. 

The  most  pronounced  case  was  that  of  the  C„Fn , - i-0oHn  0 mixture 

f 16  — o lo 

containing  12.13  volume  per  cent  fluorocarbon.  For  this  mixture  at' room 
temperature  the  liquid  rose  to  a level  approximately  8 cm.  above  the 
meniscus  of  the  bulk  of  the  liquid.  Creeping  was  observed  to  a lesser 
degree  for  mixtures  containing  up  to  UO  per  cent  C^F^,  above  which 
composition  the  phenomenon  did  not  occur. 

For  system  eight,  a vertical  creep  of  2.5  cm.  was  observed  at 
room  temperature  for  the  11.18  volume  per  cent  CgF.^0  mixture  although 
none  was  observed  in  the  mixtures  composed  predominantly  of  fluoro- 
carbon. System  seven  showed  this  phenomenon  to  a lesser  degree  in  the 
mixtures  containing  less  than  kO  per  cent  of  the  chlorofluorocarbon. 

The  systems  which  creep  up  the  glass  walls  all  show  minimums 
in  the  surface  tension.  An  explanation  based  upon  this  fact  alone  is 
not  adequate,  however,  as  the  CgF.^0  - n-C^H-ji  system  (Figure  17)  shows 
a much  deeper  minimum  than  the  C^Cl^  - i-CgH-^  system  (Figure  15)  and 
the  former  system  does  not  exhibit  this  unusual  behavior  while  the 


latter  one  does. 


V.  CALCULATIONS  FROM  VISCOSITIES 
OF  BINARY  SYSTEMS 


Free  Energy  of  Activation  for  Viscous  Flow.  Eyring  et  al.  [32] 
have  applied  the  theory  of  absolute  reaction  rates  to  the  viscosities 
of  liquids.  According  to  this  model,  flow  is  pictured  as  a molecule 
reaching  an  activated  state  as  it  squeezes  past  its  neighbors  and  then 
dropping  into  the  next  available  equilibrium  position.  In  this  theory, 
the  absolute  viscosity  is  expressed  as 


^ = distance  between  adjacent  moving  layers  of  molecules, 
a = distance  through  which  the  shearing  force  acts, 
h = Planck's  constant, 

N = Avogadro's  number, 

V = molar  volume, 

AG*  = free  energy  of  activation, 

R = gas  constant,  and 
T = absolute  temperature. 

For  ideal  solutions  in  the  thermodynamic  sense  it  has  been  suggested 
[32,  36]  that  the  viscosity  be  given  by 


(1) 


where: 


111  \ = Xa  ln  \ + h ln  V 


(2) 


h7 


t 


1*8 

Since  the  volume  of  a solution  varies  with  the  composition,  it  seems 
more  reasonable  to  replace  equation  2 by 

ln  % V3  ■ *«  1,1  \ \ * h ^ \ Vb  • ^ 

If  equation  3 is  used  as  the  definition  of  ideal  solution  behavior 
in  the  viscosity  sense  and  (1/a)  is  set  equal  to  unity,  the  free  energy 
of  activation  for  an  ideal  solution  is  linear  in  mole  fraction,  i.e., 

For  non-ideal  solutions,  the  observed  free  energy  of  activation  will 
differ  from  the  ideal  value  by  a term  g(x)  which  is  a function  of 
concentration: 

*0*  - AO*1  - g(x).  (5) 

By  assuming  (1/a)  is  unity  Ag* , values  have  been  calculated 

s 

by  equation  1 and  the  points  in  Figures  1 through  10.  The  AG*  values 
calculated  in  the  same  way  for  the  pure  components  were  used  in  equa- 
tion h to  give  the  ideal  solution  values  Ag*1.  The  difference  between 
these  values,  g(x),  is  plotted  for  the  25°C  and  1*5°C  isotherms  for 
systems  one  through  ten  in  Figures  19  through  28,  respectively. 

It  should  be  recalled  that  systems  one  through  ten  represent 
a steady  departure  from  ideality  in  the  thermodynamic  sense.  The  vis- 
cosity of  system  one  is  essentially  ideal  and  g(x)  has  very  small  posi- 
tive values  (Figure  19)  corresponding  to  the  small  negative  volume 
changes  on  mixing.  System  two,  which  is  also  ideal  thermodynamically, 
is  far  from  ideal  in  viscosities.  Figure  20  shows  that  g(x)  has  rather 
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Figure  19. —Excess  Free  Energies  of  Activation 
for  the  C?F16  - CqFi60  System 
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Figure  21.— Excess  Free  Energies  of  Activation 
for  the  C^Cl^F,,  - CgF^O  System 


Figure  22. — Excess  Free  Energies  of  Activation 
for  the  C^C12F6  - CgF^O  System 


51 


G 


.figure  23. — Excess  Free  Energies  of  Activation 
for  the  C^Cl^F^  - C„F^  System 


Figure  2k. — Excess  Free  Energies  of  Activation 
for  the  n-C6H^  - C^Cl^  System 


■ igure  23.  Excess  i*ree  Energies  of  Activation 
for  the  C^Cl^F^  - i-CgH^g  System 
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Figure  26. — Excess  Free  Energies  of  Activation 
for  the  i-CgH^g  - cgF-j_6°  System 


Figure  27. — Excess  Free  Energies  of  Activation 
for  the  n-C^H.^  - C8Fi6°  System 
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Figure  28. — Excess  Free  Energies  of  Activation 
for  the  i-CgH^g  - System 
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large  negative  values  over  whole  composition  range.  System  five, 
on  the  other  hand,  is  not  ideal  thermodynamically  but  g(x)  has  rela- 
tively small  negative  values  (Figure  23).  The  fourth  combination  is 
illustrated  by  system  three,  which  is  far  from  ideal  both  in  its  vis- 
cosity behavior  (Figure  21)  and  in  its  thermodynamic  behavior. 

System  four  (Figure  22)  is  particularly  interesting  in  that 
it  displays  large  positive  values  for  g(x)  when  all  the  other  systems 
from  two  through  seven  have  negative  values  for  g(x) . Each  one  of  these 
six  systems  shows  expansion  on  mixing. 

It  is  immediately  apparent  that  there  is  no  simple  correlation 
between  non-ideal  behavior  in  the  thermodynamic  sense  and  non-ideal 
behavior  in  the  viscosities  of  solutions.  At  least  two  general  effects 
may  be  singled  out  from  these  studies.  The  more  obvious  one  is  illus- 
trated by  systems  six  and  seven  (Figures  2 Ij.  and  23),  where  expansions  on 
mixing  are  'accompanied  by  decreases  in  the  free  energy  of  activation  for 
viscous  flow. 

To  account  for  the  decrease  in  AG*  due  to  expansion  on  mixing, 

s 

g(x)/RT  has  been  plotted  versus  the  total  volume  change  on  mixing  in 
Figure  2 9 for  systems  six,  seven,  and  ten.  For  the  first  two  of  these 
systems.  Figures  29A  and  29B,  the  lowering  of  viscosities  of  the  solutions 
may  be  attributed  directly  and  entirely  to  the  expansion  on  mixing. 

From  the  25°C  isotherm  of  Figure  29C,  it  is  obvious  that  a second 
effect  in  the  viscosities  exists.  This  effect  corresponds  to  the  inflec- 
tions in  g (x ) for  systems  eight,  nine,  and  ten,  Figures  26  through  28, 
respectively.  In  these  three  figures  the  curves  near  the  edges  of  the 
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Figure  29. — Excess  Free  Energy  as  a Function 
of  Volume  Change  on  Mixing 
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figures,  i.e.,  the  dilute  solutions,  resemble  Figures  21;  and  25  in  which 
the  viscosity  lowering  is  attributed  to  the  volume  change  on  mixing. 

Near  the  center  of  Figures  2c  borough  28,  a second  effect  is  obviously 
superimposed  on  the  expansion  effect. 

This  second  effect  will  be  interpreted  in  terms  of  an  associa- 
tion or  local  clustering  within  the  solution.  This  effect  is  notice- 
able for  system  ten  at  1|5°C,  some  21C°  above  the  maximum  unmixing  tem- 
perature. At  lower  temperatures  this  effect  predominates  over  the  expan- 
sion on  mixing  effect  and  the  Ac-''  becomes  greater  than  the  ideal  Ag"^. 

s s 

From  Tables  10  through  12,  it  can  be  seen  that  the  volume  changes  on  mix- 
ing show  no  inflections  corresponding  to  these  increasing  values  of  g(x) . 
Similarly,  Figure  29C  shows  that  both  the  expansion  effect  (dashed  line), 
and  the  clustering  effect  (pronounced  rise  in  g(x)  at  25°C)  are  present. 

The  positive  g(x)  values  observed  over  the  entire  composition 
range  for  system  four,  Figure  22,  suggest  that  the  second  effect  is 
overwhelmingly  large  compared  to  the  expansion  on  mixing  effect.  The 
fact  that  g(x)  in  system  five  has  such  small  negative  values,  out  of 
proportion  to  the  expansion,  no  doubt  arises  from  the  interplay  of  both 
effects  in  this  solution.  Systems  four  and  five  are  similar  to  systems 
one,  two,  three,  six,  and  seven  in  that  the  thermodynamic  free  energies 
of  mixing  are  not  sufficiently  large  to  produce  partial  miscibility. 
Inasmuch  as  they  show  a dominating  value  for  the  second  effect,  these 
solutions  are  similar  to  those  solutions  (eight,  nine,  and  ten)  which 
have  large  thermodynamic  free  energies  of  mixing.  The  conclusion  to  be 
drawn  from  these  similarities  and  differences  is  that  the  second  effect 
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is  associated  with  the  thermodynamic  entropy  of  mixing  for  the  solu- 
tions. We  will  see  that  the  positive  part  of  g(x)  may  be  accounted  for 
by  allowing  ^/a  to  be  a variable  in  equation  1. 

The  Factor  1/a  as  a Variable.  By  differentiating  In  T|V  in  equa- 
tion 1 with  respect  to  1/T  at  constant  composition  one  obtains 

dln(T)V)  _ 1 5(4gVT)  , 0 3ln(X/a)  • 

o(l/T)  R "'o(l/T)  2 ”d(l/T)  ’ ^ 

Rearranging  gives  an  expression  for  the  enthalpy  of  activation  Ah’', 


3(Ag"Vt) 

■^T7tT 


ain(T)V) 
d( 1/T ) 


on  5ln(X/a) 

”<H1/t)  ' 


(7) 


From  the  experimental  viscosities,  it  was  found  that  In  TlV  versus  1/T 
gave  essentially  a straight  line  over  the  temperature  range  studied. 

This  was  true  for  the  pure  materials  and  for  the  mixtures  which  showed 
no  inflection  point  in  g(x) . Thus,  for  these  cases,  at  each  composition 
the  quantity  Bln  TlV/3(l/T)  is  a constant  independent  of  temperature. 

In  systems  eight,  nine,  and  ten,  however,  where  g(x)  shows  maxima  and 
minima  with  respect  to  change  in  composition,  Bln  T]V/3(l/T)  varies  with 
temperature  in  a manner  illustrated  by  Figure  30  for  some  data  of  system 
ten.  The  maximum  unmixing  temperature  of  23.65°C  for  this  system  occurs 
at  61.1;  mole  per  cent  isooctane  (51;  volume  per  cent  at  25°C) . As  the 
temperature  for  this  composition  is  lowered,  Bln  T|v/3(l/T)  shows  a 
rapid  rise  beginning  in  the  neighborhood  of  27°C.  At  27°C,  RBln  T|v/3(1/t) 
is  1;000  calories  per  mole  and  rises  to  87,000  calories  per  mole  at  23.7°C. 
In  the  dilute  solutions,  even  for  these  systems  (eight,  nine,  and  ten), 

Bln  T]V/3(l/T)  is  a constant  independent  of  temperature. 
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Figure  30.— R 3 in  p./d(l/T)  as  a Function  of  Temperature 
for  the  i-CgHl8  - CyF-^  System 
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As  an  hypothesis  to  explain  these  observations,  we  will  assume 
the  following  schemes  for  the  flow  process. 

I.  In  pure  liquids  or  in  solutions  where  no  exaltation  of  g(x) 
is  observed, 

- A7'  + hole  ^ (I ) 

II.  In  liquids  where  g(x)  shows  the  exaltation, 

A^  - nA*  + n holes  -*  A^  (II) 

In  I,  A^  represents  a molecule  in  position  1,  k"  represents  the  mole- 
cule in  the  activated  state,  and  A ^ is  the  molecule  displaced  to  posi- 
tion 2.  In  IIj  A^  represents  n molecules  in  some  group  or  association 
(e.g.,  in  the  incipient  two-phase  region  at  temperatures  just  above  the 
unmixing  temperature).  A^  represents  the  molecules  displaced  to  new 
positions  and  included  in  a new  group. 

The  enthalpies  of  activation  per  mole,  Ah  , in  a given  liquid  for 
process  I and  for  process  II  are  essentially  identical,  since  the  holes 
required  possess  the  same  properties  in  both  processes.  Thus,  Ah"’  will 
be  a function  only  of  composition  in  either  process  I or  II . In  the 
temperature  region  where  the  viscosity  or  g(x)  versus  composition  show 
the  maxima,  under  these  hypotheses,  the  excessively  high  temperature 
coefficient  of  viscosity,  and  of  g(x)  above  those  values  expected  in 
the  absence  of  the  grouping  effect,  must  be  exactly  compensated  by  the 
term  -2R  9ln( Va)/d(l/T)  in  equation  7.  It  is  shown  in  the  following 
paragraphs  that  this  is  true.  This  hypothesis  thus  attributes  a posi- 
tive portion  of  the  total  g(x)  to  a change  in  cluster  size  with 
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temperature.  The  shearing  force  acts  on  an  area  X which  is  larger 
when  clusters  are  present  than  when  they  are  absent. 

Classification  of  Viscosity  of  Solutions.  A more  explicit  form 
of  equation  5 is  needed  which  includes  values  of  X/a  different  from  unity. 
From  equation  1 we  may  write  for  the  solution 

ln(71  V/hN)  = AgYrt  + ln(-)2  . (8) 

So  S 3. 

If  we  then  write, 

ln(TlsVs/hN)  = xa[^°/RT  + ln(X/a)2]  + xJ^*°/RT  + ln(X/a)2] 

+ gv(x)/RT  + In  f (9) 

the  term  In  f is  defined  by, 

In  f = ln(X/a)a  - x&ln(X/a)a  - xbln(X/a)2  . (10) 

Then  the  g(x)  of  equation  3 becomes 

g(x)  = gy(x)  + RT  In  f,  (11) 

in  which  g (x)  is  that  part  of  g(x)  which  arises  from  the  change  in 
volume  on  mixing.  Equations  7,  9,  10,  and  11  serve  as  the  basis  for 
classifying  the  viscosity  behavior  of  solutions: 

1.  For  ideal  solutions  both  g (x)  and  In  f are  zero.  Thus 

g(x)  is  zero  and  the  ideal  solution  obeys  equation  3.  The  two  fluoro- 
carbons, C-^F^  CgF-^0,  which  constitute  system  one,  are  very  close 

to  this  behavior. 

2.  For  this  class  of  solution  g(x)  is  not  zero  but 

[ dln(11V)/d(l/T)  ] is  independent  of  temperature  and  (X/a)  is  independent 

X 
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of  temperature  for  any  given  composition.  Subclasses  would  be  those  for 
which  g(x)  is  positive  or  for  which  g(x)  is  negative. 

3.  The  third  type  of  mixture  is  that  in  which  [ dln(T|V)/d(l/T)  ] 
is  a function  of  temperature.  Systems  eight,  nine,  and  ten  fall  in 
this  class. 

Treatment  of  Class  3 Solutions.  At  sufficiently  high  tempera- 
tures Class  3 solutions  behave  like  Class  2 solutions  and  for  each 
->£ 

composition  AH  is  a constant  by  equation  7 since  [ dln( Va)/d(l/T) ] is 

then  equal  to  zero.  These  values  of  AH  agreed  within  a few  per  cent  of 

s 

those  calculated  by  assuming  AH"  was  ideal,  i.e., 

s 

- xa  < * *b  *b  <12> 

These  ideal  values  of  AH"'1'  were  used  in  these  calculations.  Using 
these  values  in  equation  7,  [2  dln(X/a)  /3(1/T)]  may  be  evaluated. 
According  to  equation  10, 


3lnf\  = /S ln(X/a)s\ 

^4'  Ndw  L- 


(13) 


since  b]_n(k/a)/S(l/T)  is  assumed  to  be  zero  for  the  pure  compounds. 
This  last  is  taken  to  be  true  since,  as  with  Class  2 solutions, 
6ln(flV)/3(l/T)  is  independent  of  temperature  for  the  pure  components. 

Substitution  of  equation  13  into  equation  7 and  integrating 
between  T^  and  T^  at  a given  composition  gives, 


1 *2 

in  -rr- 


ln- 


m. 


AT 

D 


--1 

T2  Tl-1. 


(14) 
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Since  the  composition  is  fixed 

m2  (Vp)2  ^ 

"(,,!v)1  " Wp  )1  - 


(15) 


where  p*  is  the  kinematic  viscosity.  Using  kinematic  viscosities  in 

equation  lit  and  the  ideal  AH  ' values  discussed  above  gives  the  ratio 
f2 

of  shape  factors  -g-  at  two  temperatures. 

1 

The  difference  in  g(x)/RT  for  the  two  temperatures  is  given  by- 
equation  11  as 

g(x)2/RT2-  g(x)1/RT1=  gv(x)2/RT2-  g^x^RT-j.  + In  fg/f^  (16) 


It  is  shown  in  Figures  29A  and  29B  that  gv(x)/RT  at  two  temperatures, 
25°C  and  1|5°C,  are  essentially  equal  if  the  volume  change  on  mixing  is 
not  very  different  at  the  two  temperatures.  Setting  gv(x)^/RT^  = 
gv(x)2/RT2  and  combining  equations  11,  li|.,  and  15  gives 

g(x)1/RT1  = g(x)2/RT2  - In  + ( AH^/R)  [l/T2  - 1/T^  . (17) 

Referring  now  to  Figure  31  for  system  ten  and  taking  T^  = 1;50C  and 
T2  equal  to  the  various  temperatures  between  h$°C  and  the  critical  , 
unmixing  temperature,  g(x)/RT  at  U5°C  has  been  computed  from  g(x)/RT 
at  lower  temperatures  by  means  of  equation  17.  The  points  shown  on 
the  figure  are  the  values  so  computed.  It  is  obvious  that  the  theory 
accounts  very  well  for  the  difference  in  g(x)/RT  between  the  lower  tem- 
peratures and  li5°C  since  the  calculated  points  coincide  almost  exactly 
with  the  experimental  curve  at  It5°C  obtained  from  the  viscosities  alone 
at  U5°C.  The  dashed  curve  drawn  on  Figure  31  is  g^(x)/RT  at  li5°C 
estimated  from  the  volume  changes  on  mixing  at  each  composition  by  means 
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of  the  dashed  line  drawn  on  Figure  29C.  The  same  computation  has  also 

been  made  for  the  25°C-curve  of  system  eight  (Figure  32).  The  points 

calculated  from  the  data  at  25>°C  fall  on  the  curve  for  b%°C. 

This  agreement  establishes  the  idea  that  the  exaltation  in  the 

viscosity  of  these  mixtures  arises  from  a variation  in  (A/a)  with  tem- 

s 

perature  and  with  composition.  The  small  magnitude  of  this  variation 
is  illustrated  by  values  in  Table  13.  An  increase  of  only  about  6 per 
cent  is  sufficient  to  account  for  the  greatest  increase  in  viscosity  at 
25°G  above  the  value  expected  from  the  volume  change  effect. 

A value  for  the  viscosity  at  the  critical  unmixing  point 
(T  = 23.65°C,  X = 0.6lli  isooctane)  for  system  ten  has  been  estimated 
by  extrapolating  the  value  shown  in  Table  lU  of  In  f^/f^o  at  this  com- 
position. The  kinematic  viscosity  at  the  critical  point  is  0.6U5  centi- 
stokes.  A/a  at  23.63°C  and  O.61I4  mole  fraction  isooctane  is  l£  per  cent 
larger  than  its  value  at  h5°C.  From  the  dashed  curve  on  Figure  31 

a value  for  f/f  at  1|5°C  may  be  estimated  as  1.0U  at  0.611;  mole  fraction 
o 

isooctane,  so  that  f/f  at  the  critical  unmixing  point  is  approximately 

1.19. 

Vapor-liquid  equilibria  of  n-C^  F16  with  isooctane  have  been 
reported  [69]  at  the  three  temperatures,  30°,  £0°,  and  70°C.  The  total 
excess  free  energy  and  entropy  of  mixing  at  these  temperatures  are  posi- 
tive. At  the  lower  temperatures  the  excess  entropy  as  a function  of 
composition  resembles  the  negative  of  the  curve  for  g(x)  at  1|5°C  on 
Figure  28.  The  increase  in  A/a  with  decrease  in  temperature  is  accom- 
panied by  an  increase  in  order  among  the  molecules  in  the  solution. 


g(x)/RT 


Figure  31. — Excess  Free  Energy  Function 


for  the  i-CgH^g  - C 


3^F^6  System 


g(x)/RT 
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Figure  32. — Excess  Free  Energy  Function 

for  the  i-CgH^  - CgF.^0  System 
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TABLE  13 

Shape  Factors  For  the  i-CgH^g  - CqF. 

and  i-C0H„  n - C_F. ✓ Systems 
— o ±o  ( lb 

Mole  Fraction 

2 ln[f25/fu5] 

C8F16° 

System  8:  CgF^O  with  i-^gH-^g 

25 

+0.0252 

36 

+0.0444 

50 

+0.0306 

64 

+0.0057 

75 

-o.oo5i 

--C8Hl8 

System  10:  i-CgH-^g  with  C^F^g 

45 

+0.0416 

53.1 

0.0790 

61 

0.1066 

70 

0.0660 

80 

0.0137 

f25/f45 


1.013 

1.022 

1.015 

1.003 

0.99 


1.021 

1.040 

1.055 

1.034 

1.007 
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TABLE  ll| 

Shape  Factors  at  the  Critical  Composition 
for  the  i-CgH^g  - C^F.^.  System 


T,  °C 

2 ln[fT/f^o] 

fT/,fli5° 

System  10  at  the  Critical  Composition, 
61. ii  Mole  Per  Cent  Isooctane 

23.65 

(+0. 26ll|) 

1.15 

23.7 

0.2459 

1.130 

2l| 

0.1777 

1.093 

25 

0.1066 

1.055 

30 

0.0li6l 

1.023 

U5 

0.0 

1.00 
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This  conclusion  is  also  obvious  from  the  behavior  of  the  heat  capacity 
of  this  solution  in  the  neighborhood  of  the  critical  unmixing  tempera- 
ture [52] . 


1.03  and  occurs  at  approximately  75  mole  per  cent  methanol. 

Treatment  of  Class  2 Solutions.  These  mixtures  are  those  for 
which  the  second  viscosity  effect  (clustering)  is  not  a function  of 
temperature.  Systems  two,  three,  six,  and  seven  fall  in  the  subclass 
for  which  g(x)  is  negative  for  a positive  volume  change  on  mixing. 
Systems  eight  through  ten  would  also  fall  into  this  subclass  at  suffi- 
ciently high  temperatures  or  in  the  dilute  solution  regions.  The  sub- 
class, for  which  g(x)  is  greater  than  would  be  expected  from  the  expan- 
sion on  mixing,  is  represented  by  systems  four  and  five.  Both  positive 
and  small  negative  g(x)  value  may  occur  in  the  latter  subclass. 

In  mixtures  such  as  system  seven,  for  which  g(x)  is  not  a func- 
tion of  temperature  (Figure  25) 


Data  [37]  for  methanol-toluene  solutions  fall  into  this  class 
and  fit  equation  17.  The  maximum  value  for  f^/f^  for  this  system  is 


(18) 


and  from  equation  11 


x 


= -Rlnf  = -s  (x) . 

vv 


(19) 


Substitution  in  equation  11  then  gives 


g(x)  = gv(x)  + T sy(x)  = hy(x) 


(20) 
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so  that  even  though  In  f is  not  necessarily  zero,  g(x)  may  depend  upon 
the  volume  change  for  such  systems. 

System  four  has  small  positive  volume  changes  accompanying  mix- 
ing. Nevertheless,  g(x)  for  this  system  is  positive  at  all  compositions 
(Figure  22) . The  positive  portion  of  g(x)  is  far  greater  than  the  nega- 
tive  value'  of  gy(x)  in  this  system.  If  it  be  true  in  this  case,  as  was 
suggested  above  for  system  ten,  that  large  values  of  ^/a  (and  consequently 
large  values  of  f by  equation  10)  mean  abnormally  low  entropies  of  mixing, 
it  should  be  expected  that  system  four  would  exhibit  negative  excess 
entropies  of  mixing.  The  thermodynamics  of  mixtures  of  CgF^O  with 
C^Cl^F^  should  be  interesting  in  this  respect. 

It  has  been  postulated  throughout  that  the  enthalpy  of  activa- 
tion is  the  same  for  a given  mixture  whether  the  process  be  I or  II. 

By  equation  7 this  hypothesis  leads  to 

[S  ln(Va)  /d(l/T)]  = C (21) 

o a.  A 

as  a characteristic  of  Class  2 solutions.  C^.  is  a constant  independ- 
ent of  temperature  but  a function  of  composition.  In  the  treatment  of 
Classes  2 and  3 we  have  taken  C^.  arbitrarily  as  zero,  whenever 
[d  In  7|V/d(l/T)]  is  independent  of  T.  If  C is  not  zero,  then  the 
base  Ah'  appearing  in  equation  (ll|)  is  actually  AH  +2  RCx.  The  treat- 
ment of  Class  3 on  this  basis  would  not  be  altered. 


VI.  EXAMINATION  OF  THEORIES  OF  LIQUID  VISCOSITY 


In  examining  the  various  theories  of  viscosity,  use  has  been 
made  of  52,  molecules.  These  molecules  are  listed  in  Table  15  together 
with  the  major  sources  of  experimental  data  used. 

Temperature  Dependence  of  Viscosity.  Of  the  numerous  equations 
that  have  been  advanced  to  describe  the  viscosity  behavior  of  liquids, 
the  one  originally  proposed  by  Reynolds  [8l]  in  1886, 

T\  = a eb//T  (22) 

is  the  most  frequently  used.  Eyring  et  al.  [32]  have  attempted  to 
interpret  the  parameters  a and  b in  terms  of  the  "absolute  reaction 
rate"  theory  of  kinetics.  The  Eyring  theory  has  been  popular  simply 
because  it  gives  a fairly  good  empirical  fit  to  the  experimental  data. 
The  agreement,  however,  is  far  short  of  the  precision  available  in 
viscosity  measurements.  When  In  7]  versus  i is  plotted  over  a wide  tem- 
perature range,  curves  are  obtained  for  practically  all  liquids.  These 
curves  have  increasing  slopes  as  the  temperature  decreases.  In  one  of 
the  attempts  to  obtain  better  agreement  with  experiment,  Jobling  and 
Lawrence  [1+9  j 50]  have  taken  a more  exact  form  of  Eyring' s theory  which 

contains  the  volume  dependence  as  well  as  that  of  Andrade  [2]  and 
_ 1 

proposed  that  In  H versus  ^ is  linear  only  at  constant  volume.  This 
frequently  used  idea  [22,  65] , known  as  the  free  volume  approach, 
results  in  two  difficulties.  The  first  difficulty  is  that  this  approach 


70 


71 


TABLE  15 

Sources  of  Data  for  52  Liquids  Used  to  Test 
- Viscosity  Theories 

Density 
references 

Simple  Molecules 


Ar 

12, 

8, 

62, 

60,  58,  56 

82, 

23, 

21 

Kr 

12 

21, 

61 

Xe 

56 

21 

Hg 

1*5, 

53, 

71 

72, 

h6 

N2 

8, 

62, 

60 

°2 

12, 

62, 

8 

F2 

78 

78 

CO 

62 

cs2 

26, 

b6 

29, 

25 

CHC13 

25, 

70 

92, 

29, 

70 

CC1U 

26, 

k6, 

62, 

25,  h7,  70,  33 

29, 

25, 

70,  33,  96 

CBr, 

4 

30 

30 

Paraffins 

1, 

12, 

8, 

16,  86,  62 

23, 

1 

°2H6 

1, 

86, 

1 

C3H8 

1, 

86, 

85 

1 

• 

2-cuhio 

1, 

86 

1 

i-cuHio 

1 

1 

w i , Viscosity 

Molecule  . J 

references 
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TABLE  15  (Continued) 


Molecule 

Viscosity 

references 

Density 

references 

Paraffins  (Continued) 

2-°5Hi2 

1,  hi 

i,  55 

i-c5Hi2 

V 

1,  U6,  55 

neo-C^H12 

1,  89,  73 

1 

1,  ii7 

1,  55 

--C7H16 

1 

1 

— -C8H18 

1 

1 

2 Me  heptane 

31 

31 

3 Me  heptane 

31 

31 

1+  Me  heptane 

31 

31 

3 Et  hexane 

31 

31  • 

2,2  di  Me  hexane 

31 

31  , 

3,3  di  Me  hexane 

31 

31 

2,U  di  Me  hexane 

31 

31 

2,5  di  Me  hexane 

31 

31 

3 Me,3Et  pentane 

31 

31 

2,2,3  tri  Me  pentane  31 

31 

2,3,3  tri  Me  pentane  31 

31 

2,3,1j.  tri  Me  pentane  31 

31 

n-C_Hno 

1 

1 

- 9 20 

--C10H22 

1,  70 

1,  70 
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TABLE  15  (Continued) 

Molecule 

Viscosity 

references 

Density 

references 

Cyclo  paraffins 

cy-0sHio 

1 

1 

cy-C6H12 

33,  b7 

95,  33,  27 

Me  cyclohexane 

1 

' 27 

Aromatics 

benzene 

1,  38,  62,  60,  U7,  33 

33,  1 

toluene 

1 

1,  29 

Fluorocarbons 

2-c5Fi2 

78 

78 

-■C6Flii 

78 

78 

-_C?Fl6 

80 

GO 

•\ 

o 

CO 

cyclo-CsF10 

78 

Other  Molecules 

• 

h2o 

1±5,  53,  h6 

JU6 

Methanol 

b6 

29 

Ethanol 

92 

29 

acetone 

Ii6 

29,  ii6 

ii  Ethylether 

25,  U6 

25 

chloro  benzene 

25 

25 

7b 


is  based  upon  the  equation  for  viscosity, 

T|V1/3  - a eo/VT 


(23) 


in  which  V is  the  molar  volume.  That  this  is  not  the  correct  relation- 
ship is  easily  verified  by  plotting  In  (Tl/p1^)  versus  p/T  which  gives 
curves  similar  to  those  for  In  T|  versus  1/T  described  above.  The  second 
difficulty  is  that  for  at  least  one  liquid  (C02)  this  approach,  with  the 
data  of  Jaeger  [U8],  gives  negative  values  for  the  activation  energy 
for  flow. 

Following  a similar  line  of  reasoning,  Doolittle  and  Doolittle 
[21;]  have  suggested 

ln  = a + TflT T (2W 

o 

where  Vq  is  the  volume  of  the  liquid  extrapolated  to  0°K.  This  equa- 
tion, in  terms  of  statistical  treatments  of  the  hole  theory  by  Weymann 
[9b]  and  Cohen  and  Turnbull  [18],  was  incorporated  by  Macedo  and  Lito- 
vitz  [59]  into  a "hybrid  equation," 


bV  /(V-V  ) £E\  /RT 

-n  o'  ' o/  vis' 

i|  = a e e 


(25) 


where  the  last  tern  is  taken  from  Eyring's  form  of  equation  22, 


T|  = a e 


AE*  /RT 
vis' 


(26) 


While  neither  equation  2l|  nor  equation  26  with  constant  parameters  give 
good  correlations  over  extended  temperature  ranges,  equation  25  is 
some  improvement  since  it  contains  an  additional  adjustable  parameter. 
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Following  Macedo  and  Litovitz,  b is  set  equal  to  unity  in  equation  25, 

V is  replaced  by  A and  B values  [78]  from 

V"2  = A - BT  (27) 

(an  equation  satisfied  within  approximately  0.1  per  cent  in  density 
in  these  temperature  ranges  by  all  of  the  liquids  used),  and  equation  25 
becomes 


In  T]  = in  aQ  + {[1  - (B/A)T]_1//2-l}"1+  ^^/RT  . (28) 

B -1  /2  i 

Plots  of  [In 71+  fl  - (1-T-)  } J']  versus  1/T  imply  negative  activa- 

tion energies  which  vary  with  temperature  for  19  of  the  25  compounds 
tested.  If,  following  Weymann's  argument  [91;],  aQ  is  considered  to  be 
temperature  dependent  and  given  by 


_ / RT  \1//2  (2  mkT)1/2 

° = ter)  “”v2^  * 


(29) 


two  values  of  may  be  obtained.  One  is  obtained  from  the  inter- 

v2/3,  V0 


cept  and  the  other  from  the  slope  of 


In  (T) 


versus 


1/T.  The  two  values  obtained  for  n-hexane  were  +2.7  kcal./mole  and 
-9-6  kcal./mole,  respectively..  Writing  aQ  = a'T  also  gave  negative 
4e“ 


vis 


values  by  equation  25.  As  a final  alternative  to  maintain  the  fit 


achieved  by  equation  25,  ^'Y±s  is  interpreted  by  the  authors  according 
to  equation  26  for  some  compounds  as 


m,*  _ / d In  7) 

^vis  - l^TITTy 


(30) 
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and,  in  other  cases  as 


AS*  = 

vis  \sD7TTjv  * 


(31) 


These  values,  however,  can  differ  by  a factor  of  five  [26,  1*9]. 

Various  equations  have  been  suggested  by  many  authors  (e.g., 
[32,  1*5,  53,  60,  87,  88,  90,  91]),  which  may  be  included  under  the 
general  forms 


In  T|  pn 

= a + £ 

(32) 

In  T|  Tm 

= a + T 

(33) 

in  which  n and  m have  values  which  vary  widely  from  one  author  to 
another.  Since  particular  ones  of  these  equations  provide  better 
fits  than  do  others  for  various  compounds,  no  general  correlation 
will  be  found  in  these  forms  and  they  should  be  viewed  as  empirical. 

Best  fits  of  the  experimental  data  by  equations  32  and  33  were  found 
in  this  work  to  require  widely  varying  value's  of  n or  m even  within 
the  homologous  series  of  normal  alkanes. 

In  the  present  work  it  was  observed  that  3 In  p/d(l/T)  varied 
more  widely  with  temperature  than  did  3 in  Tl/S(l/T).  Since  **  = T)/p  ^ 
this  suggested  plotting  In  Tip  versus  1/T.  This  strictly  empirical  plot 
provides  a marked  improvement  in  linearity  over  that  obtained  from 
equation  26.  This  linearity  was  found  to  be  true  even  for  methanol  for 
which  Innes  [1*5]  reports  that  m = 0 in  equation  33. 

The  point  to  be  made  here  is  that  viscosity  varies  much  more 
rapidly  with  temperature  than  does  density  so  that  over  small  temperature 
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intervals  the  effect  of  density  variation  is  minor.  Simple  relations 
between  these  two  quantities  should  be  viewed  as  empirical  and  for 
data  at  atmospheric  pressure  or  low  pressures 


is  a better  two-constant  empirical  equation  than  any  in  current 
widespread  use.  Equation  3k  applies  over  the  whole  liquid  range  for 
a particular  compound. 

Corresponding  States  for  Pure  Materials.  Reed  and  McKinley  [79] 
have  shown  that  the  A and  B of  equation  27, 


which  describes  liquids  at  temperatures  below  T = 0.8T  , may  be  related 
to  E , r , m,  and  n of  a Leonard- Jones  bireciprocal  pair-potential 
energy  function, 


m = attraction  exponent, 
n = repulsion  exponent, 

E = depth  of  the  pair-potential  energy  minimum,  and 
r = distance  between  molecular  centers  at  E^r. 

-1/2 

equals  the  molar  volume  extrapolated  to  absolute  zero, 


In  T)p 


(3U) 


V~2  = A - BT 


where : 


(36) 


Also, 


(37) 
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where  k = Boltzmann's  constant  and  y and  KY  are  constant  Tor  a given 
m and  n.  If  the  following  conditions  are  fulfilled,  it  should  be  pos- 
sible to  use  A and  B to  reduce  the  viscosities  and  obtain  a universal 
function  obeyed  by  all  liquids: 

1.  The  potential  of  interaction  is  the  sum  of  pair  interactions. 

2.  Pair  interactions  may  be  written  in  the  same  form  (e.g., 
equation  35  with  fixed  values  for  m and  n) . 

3.  The  partition  function  may  be  factored  into  translational, 
vibrational,  rotational  and  internal  parts,  and  the  trans- 
lational part  is  the  only  one  that  depends  upon  the  volume 
per  molecule. 


This  suggests  that  by  means  of  dimensional  analysis  the  reduc- 
ing factors  may  be  found  for  the  various  properties  of  the  liquid. 


The  reduced  property  obtained  should  then  be  a universal  function  of 
the  reduced  temperature.  Such  calculations  have  been  quite  successful 
in  correlating  thermodynamic  and  transport  gas  phase  properties  for 
these  compounds  [79,  6],  A and  B values  were  calculated  from  liquid 
densities  and  are  given  in  Table  16.  Other  values  were  taken  from 
Reed  et  al.  [79,  6,  17] . 


Absolute  viscosity,  T],  has  units  of  mass/(length  x time)  which 
means  that  the  reducing  factor  is  [ (m  E^)1^2  r*-2]"1.  Using  A and  B 
from  equations  36  and  37,  the  reduced  absolute  viscosity,  T]  becomes 


= 71  (BA) 
r Al/3  m1/2 


(38) 


where  M is  the  molecular  weight.  ^ plotted  against  a reduced  temper- 
ature, 


(39) 
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TABLE  16 


Liquid  Density  Parameters  A and  B 
from  Equation  27 


Molecule 

A, 

10  ^(mole/cc.)^ 

B, 

> 0 

mole^ 

2 <V 
cc.  K 

-3- 

o 

12.2598 

1*7.2516 

i-Vio 

2.0512 

3.86182 

-_C10H22 

0.1*221 

0.51*05 

cyclo-C^H10 

2.0173 

3.02752 

Toluene 

1.1*371 

1.88001 

Methanol 

10.3581* 

11*.  5256 

Ethanol 

U • 7 866 

6.31363 

Diethyl ether 

1.8081 

3.0061*1 

Acetone 

3.1*029 

• 5.29369 

chci3 

2.7151* 

3.95722 

os2 

U . 6667 

6.52388 

cn 

Hg(0-70°C) 


1.2  3 
50.5116 


1.16 


16.6053 
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should  give  one  curve  if  the  above  conditions  are  fulfilled.  Actually, 
a series  of  curves  was  obtained.  The  normal  alkanes  formed  a family  of 
curves  which  were  crossed  by  the  more  symmetrical  hydrocarbon,  CCl^, 
and  the  fluorocarbon  curves.  Reduced  kinematic  viscosities, 


»r  . ,(f)1/2  , 


(Uo) 


plotted  against  from  equation  39  gave  an  almost  identical  pattern 
of  curves. 

A number  of  other  combinations  were  attempted  but  in  each  case 
essentially  the  same  pattern  of  curves  was  obtained  as  resulted  when 
was  plotted  versus  T . Some  of  these  reduced  properties  plotted 
versus  Tr  were  [V(p  )1//2]r  after  Thomaes[89],  [ 'H V ] ^ , and  [TlV2]^. 
Collins  and  Raffel's  [20]  statistical  treatment  of  transport  leads  to 

T]2  = 0.510  (pA1/3)(-2.1jTI  - kT  + pV).  (i|l) 

Vap 

Reduction  by  A and  B gives 


1]2/(TMA2//3)  = f (TB/A)  (i|2) 

where  f should  be  a universal  function.  T]2/(TMA2//3)  versus  TB/A  once 
again  gave  the  pattern  of  curves  discussed  above.  The  dimensionless 
reduced  viscosity  of  equation  1|2  is  similar  to  that  of  equation  38  with 
E replaced  by  kT. 

The  molecules  having  curves  crossing  the  majority  of  the  other 
curves  in  the  above  treatments  were  those  which  have  relatively  short 
liquid  ranges.  Plots  of  T)  , p-^,  as  well  as  T)  and  p.,  were  made  versus 


both 


T -Ti 


m 


T, -T 
b m 


and 


T -Tc 

_T'  where  the  subscripts  b,  m,  and  c denote  the 
c m 
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normal  boiling  point,  normal  melting  point,  and  the  critical  temperature 

T -T 

[lj  7],  respectively.  The  Ti^  versus  - — ^ plot  gave  the  only  set 

c“  m 

of  parallel  curves  obtained  for  the  molecules  tested.  Attempts  to  cor- 
relate the  constants  required  to  move  these  curves  together  with  proper- 
ties of  the  compounds  were  unsuccessful. 

In  most  theories  of  viscosity  the  exponent  factor  b in  equa- 
tion 22  is  considered  to  be  proportional  to  an  energy  as  shown  in  equa- 
tion  26.  This  energy,  , was  plotted  versus  the  energy  factor  A/B 

in  these  reducing  parameters.  A single  curve  resulted  for  the  normal 
paraffins  and  many  other  organic  compounds  (Figure  3h) . The  most  not- 
able departures  from  this  curve  are  those  molecules  with  the  exception- 
ally high  values.  These  molecules  are  the  fluorocarbons  and  the 

highly  symmetrical  molecules  such  as  neo-pentane,  cyclohexane,  CCl^,  and 
benzene.  These  results  suggest  that  the  mechanism  for  viscous  flow 
involves  an  additional  rotational  energy  requirement  not  found  in 
molecules  with  central  force  fields.  On  the  basis  of  collision  theory 
for  momentum  transfer,  both  translational  and  rotational  energy  of 
molecules  may  be  exchanged.  Rotational  energy  exchange  and  transla- 
tional-rotational exchange  will  be  of  significant  magnitude  only  if 
the  force  fields  of  a pair  of  interacting  molecules  (or  at  least  one  of 
the  pair)  are  acentric.  The  force  fields  of  unbranched  hydrocarbons 
are  quite  central  while  those  of  the  fluorocarbons,  CCl^,  neo-pentane, 
benzene,  etc.,  are  acentric. 
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Corresponding  States  for  Binary  Mixtures.  In  a series  of  papers 
[11,  12,  9,  10,  89],  Boon,  Thomaes  et  al. , have  measured  relative  viscos- 
ities of  simple  molecules  (Ar,  Kr,  Xe,  CH^,  CD^  and  0^)  and  suggested 
a reduced  relationship  for  such  molecules  and  their  mixtures  obtained  in 
a manner  similar  to  equation  38. 

These  authors  propose  that 


In  = a + - , 


(43) 


be  reduced  to 


where 


M-r  = n/Z 


Z = (E*  r"2/m)1/2 

and  p-  is  a universal  function  of  T , where 
r r’ 

T = kT/E*. 

r 

Substituting  44  into  1+3  and  differentiating  with  respect  to  l/T 
implies  that 

, b 

b r 

T - T * 
r 


OUi) 


(U5) 


m 


(U7) 


where  b^  is  a universal  constant.  Then  from  equations  43 > 44,  and  47, 


and 


In  = a - In  Z + |- 

r 


a - In  Z = constant  for  all  substances. 


(48) 


m 


By  equations  36,  37,  and  45 

a = In  [(BM)^2/A^^]  + constant. 


(5o) 
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A plot  of  "a"  (obtained  from  *a  values  extrapolated  to  1/T  = 0)  versus 

1/2  1/3 

(EM)  ' /k  J was  constructed  with  a line  of  unit  slope  drawn  to  give 
the  best  fit  to  the  Ar,  Kr,  Xe,  and  CH^  points  (Figure  33). 

From  equations  1*6,  1*7,  and  37 

b = E1'  = (constant)  (|) . (5l) 

A plot  of  "b"  (obtained  from  d In  M>/d(l/T))  versus  (5)  was  constructed 

D 

and  the  best  straight  line  was  drawn  from  the  origin  through  the  Ar, 

Kr,  Xe,  and  CH^  points  (Figure  3l*) . 

From  knowledge  of  A and  B,  Figures  33  and  3l*  should  give  the 
a and  b values,  respectively,  which  permit  calculation  of  the  viscosities 
for  the  simple  molecules  by  equation  1*3. 

Prigogine  [75]  has  shown  that  for  his  average  potential  model, 
the  average  values  (E  ) and  (r  ),  for  the  species  in  a binary  mixture 
are  given  by 


<E*> 


_ [XlElrf  ^ + *Krf] 


6 2 * *6' 


n/(n-6) 


2E*rf  ♦ 


V*nf/(n-6) 

+ X2E2r2  J 


(52) 


<r*> 


lElrl 


-)*11 

T‘2ji12ri2 


X:E,rl'  + 2X.X.E;„r:„  + X))E"r 


2 


l/(n-6) 


_2„* 


9*  ttc 


-XlElrl  + 2XlX2E12ri2  + X2E2r2  - 


(53) 


for  the  Lennard-Jones  (6-n)  potential.  Subscripts  1 and  2 imply 
interactions  between  pairs  of  like  molecules  and  subscript  12  implies 


-o.U 


8b 


1A 


-=t 


•o 

o 

i 


o 

CM 

I 


CM 

CM 

I 


CM 

I 


M3 

CM 

I 


Figure  33. — Reduced  Relationship  for  "a"  of  Equation  bd> 
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(A/B) 


Figure  3h. — Reduced  Relationship  for  "b"  of  Equation  h3 
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pair-interactions  between  unlike  species.  Assuming  the  usual  combining 
rules. 


Equations  52  and  53  were  written  in  terms  of  A and  B by  equations  36 


the  various  mixtures  of  simple  molecules.  Values  of  a and  b from 
Figures  33  and  3U  then  gave  P'ca-^c  for  these  mixtures  which  could  be 
compared  with  the  experimental  values. 

The  calculated  values  for  the  Ar  - Kr  system  were  in  excellent 
agreement  with  experiment,  p versus  composition  shows  a minimum  of 
approximately  0.8  millistokes  [9]  below  a linear  function  and  the  cal- 
culated values  are  within  0.05  millistokes  of  the  experimental  value  of 
1.88  millistokes. 


systems  (2.76  and  1.68  millistokes,  respectively  [11])  represent  nega- 
tive deviations  of  1.1  and  1.2  millistokes  below  linear  (J.  versus  mole 
fraction  relationships.  The  calculated  values  predicted  positive 
deviations  of  0.5  and  0.2  millistokes  above  the  linear  relationship, 
respectively. 

The  0^  - Ar  system  exhibits  only  very  slight  positive  deviations 
in  P-  from  linearity.  The  experimental  value  of  1.96  millistokes  [28] 


(5U) 


and 


(55) 


and  37  and  values  of 


were  calculated  for 
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represents  a deviation  of  only  +0.0^  millistokes  while  the  calculated 
value  predicts  a positive  deviation  of  +0.?  millistokes. 

For  the  above  systems  of  simple  molecules,  other  than  Ar  - Kr, 
the  predicted  values  are  definitely  too  high.  For  the  benzene  - CCl^ 
and  benzene  - cyclohexane  systems  (|^  = 5.967  ms.  and  7.3hh  ms.,  respec- 
tively [ 33,  19])  the  calculations  gave  values  which  were  too  low  (devia- 
tions of  -2.7  and  -k.7  millistokes  compared  to  deviations  of  — 0 . U and 
-1.6  millistokes  experimentally). 

Lennard-Jones  (m-n)  values  for  the  above  compounds  have  been 
estimated  [79]  to  be  (6-12)  for  Ar,  Kr,  and  CH^,  (6-17)  for  0^,  and 
(7-28)  for  benzene,  CCl^  and  cyclohexane.  Changing  equations  52  and  53 
to  fit  a (7-28)  potential  function  did  not  noticeably  alter  the  calcu- 
lated kinematic  viscosity  for  the  benzene  - CCl^  system.  This  indicates 
that  the  form  of  the  potential  function  is  probably  not  the  primary 
source  of  disagreement  between  the  theory  and  experiment. 

The  primary  source  of  disagreement  is  probably  to  be  found  in 
the  assumption  that  only  the  configurational  space  part  of  the  parti- 
tion function  and  not  the  momentum  space  part  is  important  in  establish- 
ing the  corresponding  states  treatment  of  transport  properties.  Another 
defect  might  also  be  in  the  use  of  the  average  potential  model  for 


mixtures. 


VII.  CONCLUSIONS 


Liquid  viscosities  can  not  be  explained,  or  indeed  even  corre- 
lated by  any  of  the  theories  thus  far  proposed.  Any  adequate  theory 
will  apparently  have  to  include  terms  involving  energies  of  rotational 
interaction  [7U,  63,  h3,  Ul]  (most  noticeable  for  the  molecules  with 
acentric  force  fields  such  as  the  fluorocarbons  and  the  quasi- spherical 
hydrocarbons).  This  suggestion  is  supported  by  the  larger  increases  in 
the  temperature  coefficient  of  viscosity  with  branching  for  the  smaller 
hydrocarbons  (e.g.,  n-C^H^,  and  neo-C^H^)  as  compared  to  the  absence 
of  this  increase  with  branching  for  the  larger  hydrocarbons  in  which 
the  chains  can  coil  (e.g.,  the  isomeric  CgH^g  compounds  in  Table  16). 
Other  evidence  lies  in  the  fact  that  CBr^,  with  its  extremely  acentric 
force  field,  exhibits  the  highest  temperature  coefficient  of  viscosity 
of  any  of  the  52  molecules  tested.  Further,  Andrade  and  Dodd  [3]  have 
found  that  for  polar  liquids  an  electrical  field  applied  in  the  direc- 
tion of  flow  has  little  effect  on  the  viscosity  whereas  a transverse 
field  produces  a marked  increase  in  viscosity.  This  last  observation 
may  be  explained  as  an  additional  rotational  interaction  when  the  mole- 
cules are  either  aligned  perpendicular  to  the  direction  of  flow  or  have 
formed  clusters  about  ions  [3].  A field  applied  in  the  direction  of 
flow  would  be  expected  to  have  little  influence  on  the  viscosity  as 
the  molecules  would  tend  to  align  in  the  direction  of  flow  and  thus 
not  increase  any  rotational  interaction.  Finally,  the  only  mixtures 
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for  -which  the  corresponding  states  treatment  predicted  correct  values 
of  the  viscosity  were  those  in  which  the  molecules,  Ar  and  Kr,  have 
truly  central  force  fields. 

The  fact  that  corresponding  states  treatments,  even  for  con- 
formal systems,  gave  results  which  deviated  further  from  experiment 
than  did  a-  simple  linear  relation  for  |Jt  in  terms  of  mole  fraction 
indicates  that  the  use  of  the  configurational  partition  function  is 
inadequate.  It  appears  that  the  momentum  integral  must  be  included  to 
account  for  rotational  and  translational-rotational  energy  exchanges. 

The  correctness  of  the  reducing  parameters,  A and  B,  used  in  this  work 
is  shown  by  the  fact  that  they  effectively  correlate  other  properties 
(including  gas  viscosities  [6])  for  these  molecules.  The  idea  that 
the  primary  difficulty  lies  in  the  partition  function,  rather  than  in 
the  potential  used,  is  further  supported  by  the  fact  that  Lennard- Jones 
(6-12)  and  (7-28)  potentials  gave  essentially  identical  values  for 
p-  , . Similarly,  a Sutherland  potential  applied  [39]  to  the  theory  of 

Born  and  Green  [13]  gave  results  which  were  at  least  as  good  as  a 
Lennard-Jones  (6-12)  potential.  It  is  obvious  that  a Sutherland 
potential  contains  an  unrealistic  repulsive  potential.  While  the 
repulsive  portion  of  the  pair  potential  is  probably  the  most  important 
in  liquid  viscosities,  correlations  with  the  attractive  portion  [35>  3k, 
83]  should  not  be  overlooked. 

Structural  changes  that  accompany  temperature  changes  [6k,  35] 
in  the  pure  liquids  and  that  accompany  both  composition  and  temperature 
changes  in  mixtures  (Class  3 solutions  in  Chapter  V)  must  be  incorporated 
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in  any  new  theory.  Eyring's  extension  of  his  theory  to  include  such 
structural  changes  [68]  is  quite  inadequate  in  that  it  is  based  upon 
hypothetical  deviations  from  equation  22  which  are  opposite  in  direction 
and  of  a different  form  from  those  actually  observed  [1,  U6,  92].  Low 
angle  X-ray  studies  [ill]  are  in  agreement  with  the  structural  changes 
proposed  in  Chapter  V for  system  ten. 

Future  viscosity  studies  should  include  the  fluorocarbons 
because  comparison  of  these  compounds  with  hydrocarbons  and  other 
materials  clearly  shows  the  inadequacies  of  the  various  theories. 
Viscosities  of  symmetrical  molecules  with  extremely  short  liquid  ranges 
should  be  studied  (e.g.,  2,2,3, 3-tetramethylbutane  and  acetylene). 

The  forms  of  the  curves  for  the  surface  tensions  of  mixtures  as 
a function  of  composition  should  be  further  investigated  as  well  as  the 
interesting  observation  of  the  rising  of  certain  liquids  on  glass  walls 
discussed  in  Chapter  IV. 
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